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ABSTRACT 

We have begun a program to search for supernovae and other transients in the fields of galaxy 
clusters with the 2.3m Bok Telescope on Kitt Peak. We present our automated photometric methods 
for data reduction, efficiency characterization, and initial spectroscopy. With this program, we aim 
to ultimately identify ^25-35 cluster SN la (~10 of which will be intracluster, hostless events) and 
constrain the SN la rate associated with old, passive stellar populations. With these measurements 
we will constrain the relative contribution of hostless and hosted SN la to the metal enrichment of the 
intracluster medium. In the current work, we have identified a central excess of transient events within 
1.25r2oo in our cluster fields after statistically subtracting out the 'background' transient rate taken 
from an off-cluster CCD chip. Based on the published rate of SN la for cluster populations we estimate 
that ~20 percent of the excess cluster transients are due to cluster SN la, a comparable fraction to 
core collapse (CC) supernovae and the remaining are likely to be active galactic nuclei. Interestingly, 
we have identified three intracluster SN candidates, all of which lay beyond R > r2oo- These events, if 
truly associated with the cluster, indicate a large deficit of intracluster (IC) SN at smaller radii, and 
may be associated with the IC stars of infalling groups or indicate that the intracluster light (ICL) in 
the cluster outskirts is actively forming stars which contribute CC SN or prompt SN la. 
Subject headings: galaxies: clusters: general 



1. INTRODUCTION 

Because little can escape the deep gravitational poten- 
tial wells of galaxy clusters, they are among the only sites 
where one can confidently construct constrained mod- 
els of chemical enrichment. The intracluster medium 
(ICM) contains mos t of the baryons in clusters (e.g., 
iGonzalez et al.1 12007D and is met al enriched, (Fe/H) ~ 
0.3^0 (e.g.. iBalestra et al"1l2007l ). However, none of the 
proposed mechanisms to transport the metals produced 
by supernovae from the galaxies in which they lie to the 
ICM have proven adequate — unl ess large gas outflows 
from cluster members are invoked (iRenzinill 19971 see also 
discussion in lDupke fe White] |2000( ). We are left to con- 
clude either that galaxies somehow eject most of their 
metals, some basic part of the chemical evolution model, 
such as the SNe rate or yield, is incorrect, or an enrich- 
ment channel has been overlooked. 

One such channel that has recently been pro- 
posed is enrichment by intracluster stars (ICS 
whose SNe pollute the ICM in situ (ISasakil 1200 ll 
Zaritsky et~aL1 l2004t iDomainko et all 12004: iLin fc Mohd 
20041: ISivanandam et all 120071) . Because the ICS 
comprise ~30% of t he to tal cluster stellar budget 
([Gonzalez et al.l [20051 . [2007 •. a comparable fraction of 
the ICM metals should come from this population even 
if galactic winds are highly efficient. A quantitative at- 
tempt to determine whether the dominant source of met- 
als is the ICS or the cluster galaxies has been unable 
to distinguish between the two sources because of un- 
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certainties in cluster SN e rates and SN la progenitors 
([Sivanandam et aLlfeOOTf ). 

Despite the widespread use of SNe la's as a probe of the 
cosmic expansion history, the p rogenitor systems of SNe 
la re main unidentified (e.g., Hillcbra ndt fc Niemeverl 
12000?) . Understanding the origins and basic physics be- 
hind SNe la is essential, given that their use as standard 
candles rests on the assumption that their basic proper- 
ties are constant as a function of cosmic time and that 
their luminosities can be calibrated with a n unchang- 
ing e mpirical relation (see recent review by iFilippenkol 
2005). One interpretation of the shortfall of metals in 
the chemical evolutionary models is that there exists a 
population of SNe that we are missing. Because we do 
not understand the origin of la's, this possibility must be 
considered. 

Independent recent results indeed suggest that 
SN la come fr o m tw o different stellar popu l ations 
(iMannucci et all 120051 : IScannapieco fc Bildstenl 120051 : 
iSullivan et al.ll2006l ). One type is associated with star 
forming regions, which tend to be over luminous, have a 
rate proportional to the star formation rate, and appear 
within 1 Gyr of the formation of the stellar population. 
The other type is associated with old stellar populations, 
have a rate that is proportional to the total stellar mass of 
the system, and only appear many Gyrs after the forma- 
tion of the stellar population. A SN search that probes 
both young and old stellar populations must use the star 
formation history of the survey v olume to disentangle 
the two SNe la populations (e.g., ISullivan et al.l [2009 ) 
and measure the SN rate for each population. This nec- 
essary modeling leads to larger uncertainties in the in- 
ferred rates. 

Old, elliptical galaxies only give rise to SN la associ- 
ated with the long delay time population ([Sullivan et al.l 
|2006( ). which makes galaxy clusters the natural environ- 
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ment in which to constrain the rate of this type of SNe 
la. By pursuing a SN survey in low redshift galaxy clus- 
ters the modeling problem is averted because old stel- 
lar populations dominate. The star formation history of 
these galaxies can be approximated to high accuracy as a 
star f ormation burst at high redshift (e.g.. iHolden et aTl 
2005). Furthermore, the profusion of old ellipticals in 
clusters also makes it possible to measure the rate of 
these SNe even though the rate is quite low. 

There have been several modern SN surveys that have 
focused on t he galaxy cluster environment , beginning 
with that of iNorgaard-Nielsen et all (|1989l ). The first 
systematic surv ey, the Mount Str omlo Abell Cluster Su- 
pernova Search (|Reiss et al.lll998h . ran for 3.5 years and 
discovered 52 SNe, 25 of which were associated with the 
cluster targets. The survey focused on Abell clusters at 
0.02 < z < 0.08, and of the 25 cluster SN la, 14 were 
confidently determined to be SN la. This survey did 
not at tempt a SN la rate measurement. IGal-Yam et al.l 
(2002) searched Hubble Space Telescope archival images 
of nine galaxy clusters that had more than one imaging 
epoch for SNe. Despite finding only two probable cluster 
SN, they were able to put the first real limits on the SNe 
rate in clusters at both low and high redshift. This same 
group has undertaken a longer term program in low red- 
shift (0.06 < z < 0.19) galaxy clusters, as part of the 
Wise Observatory O ptical Transient Search (WOOTS; 
IGal-Yam et al.ll2007f) . Their results incl ude the discovery 
of two IC SN la (jGal-Yam et all 120031 ) and a measure- 
me nt of the SN la rate i n galaxy clusters based on six SN 
la (ISharon et al.ll200 7rj). Very recently. IMannucci et al.l 
(2007) presented a measurement of both the SN la and 
CC SN rate in local galaxy clusters, measuring the clus- 
ter SN la rate to be ~3 times higher than that in field 
E/SO galaxies. Work has also begun at higher redshift 
(|Sharon et al.ll2007a[ ). 

Motivated by all of the above, we began a long-term 
photometric monitoring campaign of ^60 X-ray selected 
galaxy clusters at 0.1 < z < 0.2 using the 2.3 m Bok 
telescop e on Kitt Peak with its 90Prime 1-degree field 
imager (| Williams et al.ll200l . We aim to measure the 
old-population, SN la rate, both in cluster galaxies and 
among the ICS. The data are also being used to study 
other transient phenomena. Here, we present the obser- 
vational design of the survey, transient detection tech- 
niques and efficiency calculations, initial spectroscopy 
and first science results. In § [21 we describe the clus- 
ter sample, survey design and scientific goals. In § [3] 
we discuss our automated image reduction pipeline and 
transient detection technique, followed by our calculated 
transient detection efficiency in § [4j We present our cur- 
rent transient catalog in § [5j along with a sample of tran- 
sients that do not appear to lie within galaxies, which 
are generally referred to as hostless. Our initial spec- 
troscopy is presented in § [6] We present evidence for a 
central overdensity of transient sources associated with 
our cluster sample in § [7] and discuss the nature of these 
sources in § HI We devote § [9] to a discussion of the in- 
triguing implications of the three IC SN candidates at 
R > ^200- Finally, in § [10] we outline our future plans. 
Where necessary, we adopt a flat cosmology with Ho=70 
km s^Mpc" 1 , fl m =0.3, and f2 A =0.7. 

2. SURVEY DESIGN AND GOALS 



The basic observational goal of this survey is to confi- 
dently identify SNe la within the galaxy clusters in our 
sample at 0.1 < z < 0.2, whether they lie among the 
ICS (hostless) or within cluster galaxies (hosted). Given 
our observational strategy and depth, we will eventu- 
ally correct the observed SN rate to the true rate, mea- 
sure whether there are significant differences among the 
hosted and hostless events, and ultimately improve the 
model of the metal enrichment history of the ICM. 

Various considerations play a role in the survey design, 
including the available field of view, image depth, survey 
cadence, and sample size. From the depth we achieve 
with modest exposure times (total 1600 sec, see below), 
and a desire to have control fields near the clusters, we 
decided to constrain the sample to redshifts between 0.1 
and 0.2. As described in more detail below, we roughly 
center our clusters on CCD chip 1 of 90Prime, the prime 
focus wide field imager on the 2.3m Bok telescope on 
Kitt Peak. We use the other three chips to measure the 
background transient detection rate, given that there is 
a ~500 arcsecond gap between them. The ~30' x 30' 
field of chip 1 corresponds to a radius of ~1.6 Mpc at 
z=0.1 and ^3.0 Mpc at z=0.2 and allows us to probe 
out beyond r2oo in all of our clusters. 

The peak magnitude of type la supernovae is Mb = 
— 19.25, with a small intrinsic dispersion of 0.2 — 0.3 mag 
(Phillips 1993). For the redshift range of our cluster sam- 
ple, this corresponds roughly to 19 < g < 20 mag. Type 
la SNe fade by ~ 1 mag in 15 rest frame days. By tak- 
ing four 400 second g band exposures with small dithers 
we easily find SN candidates down to g ~ 21.5 during 
good observing conditions (see § [4]) , allowing us to ob- 
serve SNe la's well before and after peak brightness. We 
will present our formal effective visibility time, the time 
during which a given survey image is sensitive to SNe, 
in future work as we determine the SNe la rate in our 
cluster sample. We present preliminary calculations in 

§o 

For an exposure time of 4x400 seconds per cluster field 
and a generous allocation of ~4 nights per month on 
the 2.3m Bok telescope (excluding the southern Arizona 
monsoon season), we can survey about 60 galaxy clus- 
ters monthly. In addition to spectroscopic followup of 
a fraction of our cluster SNe candidates, we are also ex- 
perimenting with multi-band imaging of each cluster field 
at each epoch to exploit recent progres s in single-epoch 
photo typing of supernovae candidates (jPoznanski et all 
2006), but do not describe the latter here. We also em- 
ploy other methods for statistical determinations of the 
rate of cluster transients (see § [TO]) ■ 

As will be discussed in § 18.11 we expect to find ^2-3 
cluster SN la per month of 90Prime imaging employing 
the above strategy. If the program runs for two observ- 
ing seasons, each spanning six months, we will discover 
^25-35 cluster SN la, with between ~5 and 10 of these 
being hostless, IC events. This survey, when completed, 
will have 4-6 times the number of cluster SN la than 
previous focused searches which have sought to measure 
the rate in clusters. Overcoming previous surveys' small 
number statistics is essential for pinning down the rate 
and understanding the long-delay time SN la population. 

Several secondary science goals can also be pursued 
with these data. For example, active galactic nuclei 
(AGN) and quasi-stellar objects (QSOs) constitute a 
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large fraction of our transient objects (see § \5\ and [6]). 
Identification of AGN/QSOs through their variability 
is highl y complementary to other AGN selection tech - 
niques fSaraiedini et al j|2003l feOOd iTotani et al.ll2005t ). 
and may be useful in clus ter fields where there is a 
known overdensity of AGN (Rudcr man fc Ebelingl 120051 : 
iMartini et al.l 12007). Another example is the use of our 
4x400 second i mage sequence to st udy short time scale 
variability fe.g.. iBecker et aLll2004f ). A final example is 
the study of rare, mysterious transient events in cluster 
cores such as that identified in Abell 267 which at the 
time could not be fit by any known supernova template 
(jStern et al.ll200l . 

2.1. Sample Selection 

We chose our cluster sample from among 0.1 < z < 
0.2 X-ray selected galaxy clusters visible from northern 
hemisphere telescopes with a preference for galaxy clus- 
ters with known strong lensing features (to increase our 
chances, albeit low, of discovering strongly lensed tran- 
sient phenomena). Clusters were preferably chosen to 
have a high X-ray luminosity, although not in any sys- 
tematic way. The cluster sample, along with the observ- 
ing epochs included here, is detailed in Table [T] The 
X-ray luminosities are those reported directly from the 
BAX online cluster data base, which attempts to homog- 
enize luminosities based on flux measurements and uses a 
Hq=50 km s _1 Mpc _1 and Q m =1.0 cosmology. We refer 
the reader to the BAX web site for details 4 . 

3. PHOTOMETRIC TRANSIENT SEARCH 

The observatio ns were conducted wi th the prime focus 
imager, 90Prime (| Williams et al.l l2004). on the 2.3m Bok 
telescope on Kitt Peak. The focal plane array consists of 
a mosaic of four Lockheed 4k x 4k blue sensitive CCDs 
(the peak QE is near 4500 A). With the pixel scale of 
0'.'45 per pixel, this detector array gives a field of view 
of about one square degree. The chip gaps are large, 
approximately 500 arcseconds on the sky. For our pur- 
poses, these gaps are actually an advantage because it 
places the other chips farther away from the cluster. We 
center our cluster fields on Chip 1 (cosmetically the best 
out of the four), giving coverage out to at least r 2 oo for 
all targets. The other three chips provide control fields. 
Here, we present transients from chip 3, in addition to 
chip 1, to quantify the background transient rate. A 
charge trap in chip 3 decreases the effective field of view 
by approximately one-eighth, but otherwise the chip is 
relatively free of defects. 

The images considered here were obtained between 
January and June 2006, during four observing runs (our 
pilot program), and one night in December 2006. The 
single December night was followed by one night of spec- 
troscopy, which we present in § [6] Several other runs be- 
yond the pilot program, which implemented this imag- 
ing plus spectroscopy strategy, were spoiled either due 
to weather or instrument problems, leaving only the 
one useable night in December 2006 where everything 
worked. There are five clusters (noted in Table [T]) whose 
only additional epoch of imaging were taken during the 
period of the Fall 2006 campaign when the 90Prime im- 
ager was operating with anomalously high readout noise. 

4 http://bax.ast.obs-mip.fr/ 



Since these images are of lower quality, they are not gen- 
erally included in this study, but are used only as the 
reference image for differencing of the good December 
data in these five clusters. 

3.1. Basic Image Calibration 

The image processing consists of various steps that we 
describe here. We do the initial image reductions us- 
ing a nearly automated IRAF script that relies heav- 
ily on the mscred mosaic data reduction system. This 
script corrects for cross talk, trims and subtracts the 
overscan region, applies an additional bias subtraction 
to remove structure seen in the bias exposures, flat-fields 
the data and rejects cosmic rays in each individual expo - 
sure using the LACOSMIC task (Ivan Dokkuml feOOl'l. 
Regions around saturated and nearly-saturated objects 
are heavily masked because residuals near poorly sub- 
tracted bright objects can be one of the largest sources 
of false detections in our transient detection pipeline. We 
calculate exposure weight maps using a combination of 
a normalized twilight flat field and the cosmic ray re- 
jection map. We edit each image header to include the 
refined world coordinate system solution obtained using 
the msccmatch task and the USNO-A2 catalog. The final 
astrometric solution has a typical rms scatter of ~0'.'2- 
073. 

To combine the individual exposures, we rescale and 
resampl e each onto the same astr ometric grid using SEx- 
tractor (|Bertin & Arnoutsl IT996) and SWarp 5 . A final 
median combined cluster image is made for each observ- 
ing epoch using the individual exposure weight maps and 
SWarp. We chose to median combine the frames to re- 
move fast-moving objects, although slowly moving ob- 
jects, those that move by less than the typical seeing at 
the Bok telescope during the course of a 4 x 400 image 
sequence (see Figure [l), remain. 

We photometrically calibrate the final images (both the 
combined and individual exposures) using non-saturated 
point sources in common with the Sloan D igital Sky 
Survey (SDSS; lAdelman-McCarthv et all (l200l ) 6 when 
there is overlap. Of the fifty-five clusters in this study, 
thirty have SDSS coverage. Otherwise, we bootstrap 
calibrations using the photometric zeropoints of SDSS 
clusters on clear nights. Every cluster was observed on 
at least one photometric night (loosely defined to be a 
night in which the Sloan derived zeropoint did not vary 
by more than 0.1 mag throughout the night). We use 
SExtractor to generate object catalogs of individual and 
combined images, which also provides a measure of the 
average image point spread function (PSF) across the 
field. See Figure Q] for the distribution of image FWHMs 
and Figure [5] for the distribution of relative zeropoints 
during our imaging runs. Both of these quantities play 
an important role in determining our transient detection 
efficiency ( § HJ) , and will be discussed briefly here. 

When measuring the image PSF and zeropoints by 
comparison with SDSS, we select only the high S/N ob- 
jects with SExtractor. For concreteness, an object must 
have a minimum of fifteen pixels (the DETECT_MINAREA 
SExtractor parameter) above a detection threshold of 
20 — a (detect_thresh parameter) in order to be con- 
version 2.15.7; http://terapix.iap.fr/soft/swarp 
6 http://www.sdss.org 
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sidered for the PSF measurement and to be used in 
the calculation of the image zeropoint. For the auto- 
matic PSF measurement, stars that met these criteria 
were further flagged to have a flux between 2xl0 5 and 
lxlO 6 counts. These flux limits were set empirically to 
minimize outliers due to galaxy contamination and near- 
saturation objects. A final PSF measurement (confirmed 
empirically for a subset of objects) was found by taking 
the output SExtractor parameter fwhmjmage for these 
objects and computing their mean FWHM using the 1DL 
function resistant_mean, which clips based on the me- 
dian and the median absolute deviation. As can be seen 
in Figure [I] the FWHM is not good, although this has 
been factored into our survey planning and has not kept 
us from the survey depth and transient detection effi- 
ciency desired (see § Q| . 

For the zeropoint calculation, we take the above high 
S/N SExtracted objects and WCS match with point 
sources from the SDSS archive (matches are made to 
within five 90Prime pixels - ^2.5 arcsec). The mean 
image zeropoint was calculated directly using the IDL 
function resistant_mean (typically hundreds of SDSS 
point sources and 90Prime counterparts were matched). 
On those nights determined to be photometric, a lin- 
ear fit was performed on the SDSS-derived zeropoint as 
a function of 90Prime imaging airmass to derive an at- 
mospheric extinction correction for those clusters with- 
out counterpart SDSS imaging. Using this, cluster fields 
without SDSS imaging were given a zeropoint based on 
the image airmass and a bright star catalog was created 
to be used as a photometric catalog for all other imaging 
epochs of that cluster. 

3.2. Difference imaging and transient detection 

We implement an automated pipeline to PSF-match 
and difference an individual epoch image from a refer- 
ence image, identify candidate variable objects, and post 
image cut-outs that are reviewed by a human using a 
combination of Perl and I DL scripts w hose core utilizes 
the ISIS software package (|Alardll200dh . 

In detail, the sequence is as follows. When a new im- 
age is obtained, the code identifies a reference image by 
searching for the corresponding archived image with the 
tightest PSF. The new image is registered, resampled 
onto the reference image grid, PSF matched, and sub- 
tracted. Residual objects are classified as candidate tran- 
sients. If the newly acquired image happens to be the 
best to date, it is used as the reference image and the 
second best image is subtracted. The image differencing 
is performed with the mrj_phot prog ram within the ISIS 
package and we refer the reader to lAlardl (|200dD for de- 
tails. We found the procedure works well when we split 
each 4k x 4k CCD into 16 subregions for subtraction 
and the PSF kernel convolution is allowed to vary over 
each subregion. We set our upper flux limit for detec- 
tion and image rescaling well below the saturation level 
(see §|4]). We utilize detect. csh and find.csh from within 
ISIS to find high variance peaks in the differenced im- 
ages and produce a first list of candidate variables. We 
then trim the list by removing objects near bad columns 
(< 25 pixels), near masked saturated stars (< 25 pix- 
els), nearby significant negative values in the differenced 
image (< 7 — a below a local sky value calculated with 
an annulus of 50-70 pixels; most likely due to a slight 



registration misalignment that produces dipoles in the 
difference image), and near another transient (< 25 pix- 
els). The latter are likely either residuals from the poor 
subtraction of a large object or astrometric variables. We 
quantify the likelihood of the rejection of real variables in 
§ 0J These cuts tend to reduce the number of candidates 
from several hundred to less than several tens in a given 
image. Despite our best efforts, some obviously spuri- 
ous transient candidates remain, and so human review is 
necessary. 

The principal long time scale variable candidates in 
this project are expected to be variable stars, active 
galactic nuclei (AGN), novae and supernovae. Since the 
primary goal of this work is to search for SNe in the 
clusters, we tried to screen out known AGN and vari- 
able stars by crosschecking with NED through a batch 
query. Nonetheless, previously unknown AGN and vari- 
able stars are a major source of contamination in this 
work, but will lead to exciting auxiliary science. 

4. TRANSIENT DETECTION EFFICIENCY 

To understand the recovery efficiency of transient 
events in our pipeline as a function of seeing and transient 
brightness, we add point sources to individual images 
and attempt to recover them with our transient detection 
pipeline with the unmodified reference image. These test 
point sources are randomly placed in images with a uni- 
formly distributed g magnitude between 17 < m g < 24. 

A typical efficiency run for a given image consists of 
50 separate realizations, each of which incorporates 30 
random amplitude, randomly placed point sources. The 
final efficiency calculation for that image is then based 
on the recovery statistics of 1500 fake point sources. Al- 
though we have not devoted the computer resources to 
measuring our detection efficiency in each of our survey 
images, we have sampled the parameter space (PSF and 
transparency) covered by our data. 

We present typical results as a function of these two 
parameters in Figures [3] and 2J The recovery rate is ~80- 
90% over a magnitude range of <? ~ 18 — 21.5 during clear 
sky conditions and the full range of seeing conditions. 
When we do not recover the transient it is sometimes due 
to its proximity to a masked portion of the image. Oc- 
casionally, a catastrophic image subtraction error in one 
of the sixteen subtraction subregions occurs. Normally 
this manifests itself in the form of a slight misregistra- 
tion of the subregion, which produces dipole residuals 
in the differenced image, and no viable transient candi- 
dates. Occasionally if a subregion has a very bright star 
in its sector, it will corrupt the whole subregion, leaving 
a poor subtraction. An examination of when this occurs 
did not reveal any trend towards the outer edges of the 
field of view, and so this is unlikely to be the cause of 
the central excess of transients that we find in § [7] In 
the absence of these relatively rare subtraction errors in 
a subregion, we recover over 90% of our fake transients 
with 18.5 < g < 21.5. In addition to recovering the vast 
majority of the transients, we also generally measure the 
correct magnitude. In Figured! we plot the difference be- 
tween input magnitude and that recovered as a function 
of input transient magnitude. At the faintest magnitudes 
where we detect transients (g ~ 22.5), the scatter is ~0.2 
mag, with no bias. The transient detection efficiency in 
chip 3 (the off cluster chip) is consistent with that of chip 
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1 to within the noise of the inidividual magnitude bins, 
and for the rest of this work we consider them identical. 

One additional experiment was performed to gauge the 
relative detection efficiency for transients in galaxy cores 
versus those in regions with no apparent galaxy host. We 
placed artificial SN with random g magnitudes between 
21.0 and 21.5 (the faint end out to which we reliably de- 
tect transients) at varying annuli from the BCG center 
in the same sample of clusters that we used to determine 
the detection efficiency as a function of image PSF. To 
be explicit, 50 artificial transients per bin were placed 
at annuli between and 1 PSF, 1 and 2 PSFs, 2 and 
5 PSFs, 5 and 10 PSFs, and continuing out to a max- 
imum radius of five hundred pixels from the center of 
the cluster BCG. Between 95% and 100% of the artificial 
transients were recovered, with no trend seen as a func- 
tion of radius. Those few artificial transients that were 
not found by our detection pipeline were not found due 
to proximity with a masked saturated star or other im- 
age defect, as determined by visual inspection. For this 
reason we do not distinguish detection efficiencies for ob- 
jects found in the cores of galaxies versus th ose with no 
host g alaxy. A similar result was obtained by lNeill et al.l 
(2006) for the Supernova Legacy Survey, although the 
WOOTS cluster search did find a significantly decreased 
detection efficie ncy in the cores of g alaxy hosts versus 
hostless events (|Gal-Yam et al.ll2003f ). 

We will use these calculated detection efficiencies in § [5] 
to estimate the approximate number of cluster SN la we 
should observe and to compare with the observations, al- 
though as we mentioned previously this treatment does 
not properly account for the varying maximum bright- 
ness and light curve of SN la. 

5. IDENTIFIED VARIABLE OBJECTS 

We present our transient catalogs for the cluster- 
centered chip 1 in Table [3] and the off-cluster chip 3 in 
Table |U These tables contain the position, epoch of dis- 
covery, the apparent g magnitude of the transient, the 
projected distance from the transient to the brightest 
cluster galaxy, and any additional comments. Such com- 
ments may include, for example, if the transient was de- 
tected in more than one epoch as a variable (repeat detec- 
tion; RD), if there is a correspondence to a NED source, 
if the transient is an IC, hostless SN candidate (see § 15. ip . 
or if we were able to get a spectrum of the object (see 
§ [6]). To get an idea of the magnitude distribution of 
our discovered transients, we plot magnitude histograms 
in Figure [HI The dashed histogram shows the transient 
magnitude histogram after removal of sources that were 
detected as variables at multiple observing epochs, and 
the dotted line shows the transient magnitude histogram 
after removal of all sources known to not be at the cluster 
redshift (whether found through NED or our own spec- 
troscopy). 

Of the variable objects associated with NED sources, 
most are QSOs. In the cluster-centered chip 1, 20 tran- 
sients were associated with known QSOs (9 in the off- 
cluster chip 3). In addition to these, three transients 
in the cluster-centered chip 1 are associated with clus- 
ter members that are known AGN or QSOs from SDSS. 
Two additional transients in the off-cluster chip 3 are 
SDSS AGN at the cluster redshift, but are ~5.2 and 7.7 
Mpc from the presumed cluster center. These AGN may 



be associated with infalling groups at the outer edges of 
the cluster, with implicati ons for the origin of AGNs in 
cluster environments fe.g.. iMartini et alll2007f ). We will 
discuss these cluster AGN further in § [8] Fourteen tran- 
sients in chip 1 (four in chip 3) were 'repeat detections'; 
several of these are known QSOs, but given our visibility 
time (see § 18. ip . some may be SN detected at different 
parts of their light curve. Three other transients in chip 1 
are associated with galaxies at the cluster redshift, none 
in chip 3. 

5.1. Transients without Hosts 

One of the primary goals of this survey is to identify IC 
supernovae. When a potential IC event is identified, we 
combine all of the imaging epochs in which the transient 
was not detected into a single deep image to search for a 
possible faint host galaxy. The images are combined as 
a weighted average (more appropriate for reaching faint 
limiting magnitudes than a median combination). We 
run SExtractor on these deep, combined images with a 
3 — (J detection threshold and minimum area of three 
pixels (MAG AUTO magnitudes are computed). We con- 
sider a transient to be an IC candidate if its Sextractor 
extraction radius does not overlap any object in the deep 
object catalog. 

A synopsis of our IC SN candidates is presented in Ta- 
ble[2] We list the position, magnitude, projected distance 
from the BCG (Rbcg), that projected distance in terms 
of T2oo 5 and our g-band magnitude limits (both apparent, 
fniim,g, and absolute, Mn m ^ gi assuming that the object 
is at the cluster redshift). To calculate the limiting mag- 
nitude for the faint host detection, we de fine the point 
source detection limit identically to that of lBecker et al.l 
( 2004) and the NOAO Archive definition of photometric 
depth, as 

mi = m - 2.5log(1.2Wa sky n) (1) 

where mo is the image zero point, W is the seeing of 
the image, a s k y is the sky dispersion around the mode, 
and n = 3 represents the 3 — a detection limit. Point 
source limiting magnitudes are given since dwarf galaxies 
at the cluster redshift (with our typical seeing) should be 
unresolved (e.g. dwarf galax ies in the Coma cluster have 
a typical size of ~1 arcsec: iKomivama et al.l[2002t ). To 
put the absolute g band detection limits in perspective, 
one can adop t the Virgo Cluster luminosity function as 
presented by iTrentham &: Tullvl (|2002f ). Assuming this 
luminosity function is representative of the clusters in 
our sample, then ^0.2% of the galaxy light in the cluster 
would come from dwarf galaxies below M g — —14. 

Apparently hostless transients may have a variety of 
explanations. First, these could be true IC SNe. Sec- 
ond, they may simply lie within hosts that are below the 
detection threshold of our current imaging, despite the 
detection limits described above. As the survey revis- 
its our cluster fields more and more, the constraints will 
tighten. Third, they may be flaring stars or novae in 
our galaxy, with precursors too faint given our current 
detection limits. Fourth, they may be CC supernovae 
associated with star formation in IC space, as suggested 
by recent observations of H a in a tidal tail in the clus- 
ter A3627 (|Sun et al.ll2007f ). Finally, they may be highly 
variable higher-redshift objects, such as AGN, that we 
only detect in high-state. 
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To resolve this confusion, we are making spectroscopic 
identification of all of our IC SN candidates a high prior- 
ity as the survey progresses. Due to the limited amount 
of spectroscopic time we have had so far, we were able 
only to attempt confirmation of the one IC candidate 
in Abell 516 (see § |6j) . Unfortunately, the resulting spec- 
trum is inconclusive due to poor signal to noise. Its mag- 
nitude at discovery was g=23.83, well below our nominal 
limit. We will discuss the number of IC SN candidates 
found relative to the total number of SN la expected in 
our survey up till now in § 18.11 and we speculate further 
in § |9] under the assumption that the three IC SN candi- 
dates found in chip 1 are indeed part of the cluster and 
not some foreground or background population. 

6. INITIAL FOLLOWUP SPECTROSCOPY 

Spectroscopic identification of a significant subsamplc 
of our cluster SNe candidates representing the complete 
range in brightness and environment is critical to our 
understanding of the relative populations of foreground 
SNe, cluster SN la, uncataloged AGN, and CC SNe as- 
sociated with the cluster. 

We were allocated five nights for followup spectroscopy 
during the final trimester of 2006 to implement our full 
cluster supernovae program. Unfortunately, as discussed 
in §[3j the 90Prime noise levels were greatly elevated dur- 
ing most of the observing term, except for the final imag- 
ing night when a grounding problem was fixed and the in- 
strument returned to nominal. Hence, we only have one 
night of imaging followed by one night of Blue Channel 
Spectrograph followup at the MMT in this initial post- 
pilot portion of the program. We present these data here 
as a demonstration of our ability to rapidly identify and 
confirm transient sources. 

For the spectroscopy, we use the 300 line grating with a 
central wavelength of 5800 A (corresponding to a wave- 
length range of 3200 - 8400 A) fo r our Blue Channel 
Spectrograph (jSchmidt et al.l[l989T ) observations on the 
night of December 19-20 2006. The l'.'O slit was used 
throughout the night. We use IRAF to do the stan- 
dard CCD processing, one dimensional spectrum extrac- 
tion, and wavelength and flux calibration. No attempt 
was made to remove potential host-galaxy contamina- 
tion from any spectrum. We intend to revisit this issue 
in future work, especially for the hosted events. 

A summary of our initial spectroscopic observations 
is presented in Table [5l Target selection is difficult to 
quantify. We chose targets based on their position on 
the sky at the given time during the night, with prefer- 
ence given to spectroscopically observing multiple targets 
in the same cluster field in order to minimize slew and 
acquisition times. We gave targets higher priority if they 
were near the cluster center (R < ^200) and slightly offset 
from their host galaxy (note that all four of our confirmed 
SN la's are within r2oo, but two of the four show no ap- 
parent offset from their unresolved galaxy host). The IC 
SN candidate Abell 516/ID11 (whose identification was 
unsuccessful) was also given highest priority. Otherwise, 
individual targets were selected more or less at random. 
Given our seeing (^2-2.5 arcsec), it was difficult to in- 
clude host galaxy morphology as a criteria for spectro- 
scopic selection. A seeing of 2 arcseconds corresponds 
to ^3.7-6.6 kpc throughout our cluster redshift range. 
To determine accurate, broad morphological classifica- 



tion, howeyer, re quires resolutions of ~1 kpc or better 
(jLotz et al.ll2004| ). All spectra are presented in Figures [7) 
[Si [TUl HU and [T^l Of the 50 unique transient sources 
(that is, transients that were not found to be variable in 
other epochs) discovered in our December data (Tablc[3]), 
we were able to observe 25 spectroscopically, although 
two of these transients were not confidently confirmed 
spectroscopically. We now briefly present the four broad 
categories of transients we found and highlight individ- 
ual interesting objects. We conclude this section with a 
discussion. 

6.1. Supernova la 

We confirmed four type la supernova (|Sand et al.1 
2007) in our initial spectroscopic followup (Figure [7]), 
one of which is clearly associated with the cluster in our 
sample (A1246/ID10). The other three are background. 
We match ou r spec tra to spectral templates taken from 
iNugent et all (|2002f ). No formal fit was done, but match- 
ing by eye was sufficient to identify SNe la and estimate 
the time since peak luminosity. 

6.2. Variable Stars 

Five of the transients targeted for spectroscopy proved 
to be stars (Figure [TTj) . 

6.3. Quasars 

We spectroscopically identified several previously un- 
known quasars selected from their optical variability, 
with a redshift range between z ~0.8 and 2.9, and we 
present their spectra in Figures [5] and [TU1 As expected, 
the transient position is always coincident with the center 
of the galaxy host position. 

6.4. Low-z Galaxies and AGN 

Our spectroscopic sample includes six low redshift 
(z < 0.3) galaxies (Figured]). The host galaxy for tran- 
sient Z0256/ID19 is a cluster member, and since there are 
no star formation signatures this galaxy may have hosted 
a SN la. The cluster redshift for Abell 136 (z = 0.15 7) is 
based on one galaxy redshift (|Struble fc Roodl[l991h and 
so we tentatively consider A136/ID7 to also be a cluster 
member. Since this host galaxy has [O II] in emission, 
indicating star formation, the nature of any hidden SN is 
more ambiguous since both CC and SN la occur in star 
forming regions. Both of these transients are coincident 
with the galaxy host core. Future followup spectroscopy 
of these host galaxies will allow us to subtract out the 
host contamination and possibly identify underlying clus- 
ter SNe in these systems. 

Two other sources show signatures of AGN activ- 
ity. Broad Mgll in A743/ID10 indicates that this is an 
AGN. We measure the redshift of RXJ821/ID15 to be 
at z = 0.22 from the Ca H and K absorption lines. At 
this redshift, the broad emission line near 2800A ( as- 
suming z = 0.22) is not at the correct wavelength to 
be Mgll. There is a point source ~1.7 arcseconds from 
the galaxy which is the source of the observed variabil- 
ity (we centered our slit on this point source and not 
the intervening galaxy, although contamination was in- 
evitable). We tentatively conclude that this single broad 
emission line is from a quasar at unknown redshift and 
that this is the source of the observed variability, rather 
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than the z = 0.22 foreground galaxy that is contaminat- 
ing our spectrum. Similar flares from background AGN 
superimposed on foreground 'h ost' galaxies hav e contam- 
inated SN surveys in the past (|Gal-Yamll2005( ). Finally, 
it should be noted that just because a host galaxy shows 
AGN activity does not mean that this is the source of 
the transient; a SNe is still possible in these systems. 

6.5. Unknown 

We have two spectra that we could not confidently 
identify (Figure fT2|). The spectrum of A516/ID10 cor- 
responds to a very faint, but IC SN candidate. This is a 
blue, relatively featureless spectrum. The other transient 
is offset from the host galaxy center with a spectrum that 
is also relatively blue, but has several broad features. We 
tentatively identify this spectrum as that of a SN Ib/c 
roughly eight days post explosion. The correspondence 
can b e judged relative to the overplotted template spec- 
trum (jNugent et al.1 12002] ) . However, due to significant 
host galaxy contamination and relatively low signal to 
noise, we categorize this identification as tentative. 

6.6. Discussion of spectroscopic results 

It is worth comparing our spectroscopic results with 
those obtained in similar general SNe and cluster SNe 
searches. We discuss three subjects in particular 
- the AGN frequency, the relative number of back- 
ground/foreground SNe to cluster SNe, and the field ratio 
of SN la to CC SN. Throughout this discussion, it must 
be kept in mind that we are dealing with small number 
statistics based on one night's worth of spectroscopy. 

AGN are a perennial contaminant in SNe surveys, 
some of which preferentially avoid SNe candidates in 
galaxy cores for this specific reason (e.g. Section 3 of 
iMatheson et all [2005) or have many epochs of imaging 
with which to identify intermittent AGN flaring. Given 
our poor seeing, many of our SNe candidates appear cen- 
tered on galaxies (including unresolved galaxies); indeed, 
two of the four SNe spectroscopically discovered were co- 
incident with point sources in the reference frame. This 
suggests that with our seeing conditions, we should not 
avoid SN candidates for followup which are in the nuclear 
regions of their host. It is difficult for us to compare our 
fraction of AGN to SNe with other surveys since very few 
publish all of their spectra or have much more involved 
SNe candidate identification algorithms. 

Of the four confirmed SNe la, one is in a cluster and 
three are background. While these numbers are sub- 
ject to incompleteness, very small number statistics, and 
survey field mismatch, it is worth comparing them to 
the other cluster SNe surveys. First, the Mount Stromlo 
Abell Cluster Supernova Search discover ed 48 SNe, half 
of wh ich were hosted by cluster members ([Germany et al.l 
l200l . Similarly, the WOOTS cluster SNe survey found 
that seven out of their twelve confirmed SNe were in the 
targeted clusters. While it appears that these other clus- 
ter SNe surveys found a slightly higher cluster to field 
SNe ratio, given our small sample size we are in agree- 
ment with this previous work. 

Finally, the ratio of field SNe la to CC SNe (three 
to one if the ambiguous spectrum from § 16.51 is a field 
SN Ib/c, although our consistency does not change re- 
gardless) is consiste nt with the spe c trosco pically com- 
plete field results of iGal-Yam et all (|2007t ), who found 



that four out of five of their field SN were SN la, al- 
though this is a difficult comparison to make since our 
survey is not spectroscopically complete. As in the flux- 
limited WOOTS survey, the effective volume probed for 
the more luminous SNe la is large enough that it out- 
weighs the fact that CC SN are more common per unit 
volume. 

From the above, it is clear that future work must fo- 
cus on intelligently selecting transient targets for spec- 
troscopic followup. For this purpose, we are obtaining 
multi-band, multi-epoch data for all of our transients 
as this survey continues in order to more accurately 
deter mine a phototype pri or to spectroscopic followup 
(e.g. iPoznanski et al.ll2006h . 

7. IS THERE AN OVERDENSITY OF TRANSIENTS 
ASSOCIATED WITH THE CLUSTERS? 

While spectroscopic confirmation is the most direct 
method of confirming an association between transients 
and the cluster, a statistical excess within the cluster is 
another method of confirming a general correspondence. 
This approach may not, however, enable a full distinction 
between certain classes of transient. For example, while 
variable stars will not correlate with the cluster position, 
both SNe and AGN will. 

To examine the statistical properties of the candidate 
transients, we place all of the clusters on a similar physi- 
cal scale by taking their published X-ray luminosities, as 
sho wn in Table [H and the Lx- Mim relation as found 
by iReiprich fc Bohringeri j2002) to calculate M 2 oo for 
each cluster. After converting M200 to our adopted 
ACDM cosmology, we then calculate r 2 oo using M 2 oo = 
200(4/37rr2ooPerit)j where p cr it is evaluated at the cluster 
redshift. 

We plot the transient projected density in the cluster- 
centered CCD chip 1 as a function of scaled radius, 
R/r 200 j m Figure [T3l Error bars are based on Poisson 
statistics. We measure the 'background' or 'field' pro- 
jected density of transients using the off cluster CCD 
chip 3 events (Tabled]). This background transient rate, 
0.45 ± 0.05 per A200 (^200 is the area enclosed within 
-^200)7 is plotted as a horizontal band in Figure [T3l Be- 
yond a scaled radius of 1.25, the cluster transient pro- 
jected density matches the background. 

The number of transients per A200 within R/r 200 = 1 
is 0.78 ± 0.08. The statistical significance of the excess in 
this inner region relative to the background rate is 3.7 a. 
In terms of the 107 unique candidates we have (within 
^200)7 these numbers suggest that we have roughly 45 
cluster- associated transients within i?/r 2 oo = 1. 

8. CLUSTER TRANSIENT POPULATIONS 

There are three plausible sources for the central excess 
of cluster transients: cluster SN la, cluster CC SN, and 
cluster AGN. We now discuss the relative contribution 
of each of these classes. 

8.1. Cluster SN-Ia 

To estimate the expected number of cluster SN la in 
our sample we utilize the standard SN la rate per unit 
stellar mass equation as calculated by a survey: 

N 

RSN^^^^ (2) 
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where N is the number of SN la discovered, Ati is the 
visibility time of the i-th observation, and Mqj is the 
stellar mass of the cluster being observed in observation 
i. The sum is over all individual observation epochs. By 
solving for N, given known visibility times, cluster stellar 
masses, and an adopted SN la rate from the literature we 
can approximately determine the fraction of our excess 
cluster transients that are SN la. 

For the SN la rat e, we adopt the rate measured by 
ISharon et all (|2007H ), 0.098 SNuM, for several reasons. 
First, their sample is the most analogous to the current 
sample - they measured the SN la rate in galaxy clus- 
ters at 0.06 < z < 0.19. Second, this rate is a relatively 
large one in comparison to other SN la rate measure- 
ments of E/S0 galaxy populations (although it is consis- 
tent with other measures within the uncertainties) and 
so provides a conservative upper limit on the the number 
of SN la expected in the current survey. Note that if we 
adopted the SN la rat e in local galaxy clusters found by 
Mann ucci et alj (|2007f ) - which is consistent within the 
uncertainties of Sharon et al. - our final SN la number 
would be reduced by two-thirds. For our purposes of es- 
timating the expected number of SN la at the factor of 
two level, the adopted rate is not crucial. 

To calculate the stellar masses of our galax y clusters, 
we uti lize the Lk,200 — A/200 relation found bv lLin et alJ 
(2004) and our M200 values calculated from the Lx-A^oo 
relation utili zed in § [71 Ad opting M/L* t K = 0.78 (see 
Figure 3 of iLin et alTl2"003h . we assign a stellar mass 
to each cluster. This calculation relies on aspects con- 
structed from different cluster samples and techniques 
and so should be treated cautiously. 

Finally, we need a measure of the effective visibility 
time of an observation, Ai, or the time during which a 
SN la in the cluster would be visible given the detection 
limits and transient detection efficiency of that image. 
We postpone a detailed calculation of the visibility time 
till future work where we will measure the SN la rate 
in low redshift clusters directly. For the moment, an 
accurate estimate of the visibility time will allow us to 
estimate the number of SN la we should have observed 
thus far. With this goal in mind, we use a highly sim- 
plified model for our detection efficiency based on our 
calculations in § [4] 



rj(m) 



m < 18.5 

18.5 < m< 21.5 

m > 21.5 



(3) 



Although this parameterization is only approximate, it 
has the advantage of being simple to implement. 

The final step is to determine the time that a SN la 
spends between 18.5 < g < 21.5 as a function of clus- 
ter redshift. We adopt a non-stretched peak absolute 
B-mag nitude for SN la of M B = -19.25 (Sulli van et all 
2006). In this estimate, we neglect the uncertainty in 
this value (^0.15 mag) and the correlation between light 
curve width and peak magnitude (SN la with a brighter 
peak magnitude have a longer rise and fall time around 
maximum brightness; parameteri zed by the 'stretch' re- 
lation, e.g. lPerlmutter et al.lll999l ). The latter leads to an 
upper limit because the SN la originating from old stel- 
lar populations, such as those in cl usters, tend to have a 
smaller stretch factor than typical (|Sullivan et al.ll2"006f) . 



To obtain g band light curves for SN la at a given clus- 
ter redshift, we begin with the multiepoch SN la sp ectral 
and photometric templates of iNugent et al.1 (|2002f ). Us- 
ing the B— band f i lter tr ansmission curves, we normalize 
the INugent et al.l (2002) spectra to produce the correct 
B magnitude at a given point on the light curve, redshift 
the spectra to the cluster redshift, and then synthesize 
the g magnitude using the g— band filter transmission 
curve. This method is analogous to that presented by 
ISharon et all (|2007b[ ) . and we show sample light curves at 
z = 0.1,0.2 in Figure [Ml As a check, we also considered 
host galaxy extinction by t aking E(B — V) = 0.2 (the 
l-a dispersion value used by Neill et al. 2006) and rerun- 
ning our analysis. Including this level of host extinction 
would cause ^30% fewer cluster SN la to be indentified 
by our survey, which for the purposes of the following 
discussion we neglect. In any case, the inclusion of host 
galaxy extinction into our calculation would increase the 
number of unaccounted for cluster transients. 

Putting all of these ingredients together, we calculate 
that N in Eqn[2]is ~10. Although this number is a coarse 
estimate for reasons outlined above, the uncertainty in 
this number is tens of percent rather than a factor of 
two. If we say that ~10 of the excess cluster transients 
are cluster SN la, then the other ~30-35 excess events 
must be either CC SN or cluster AGN. Importantly, the 
vast majority of the candidates are not SN la's. 

Because of the dominance of contamination, it is im- 
perative to classify the transients. Spectroscopy is clearly 
the preferred way to do this and from the results pre- 
sented so far we can test whether the above estimate is 
reasonable. If we divide the number of expected cluster 
SN la (ten) by the number of imaging nights presented in 
this study (thirteen) to conclude that our detection rate 
is ^0.8 cluster SN la per imaging night (the true detec- 
tion rate of cluster SN la is slightly higher because of 
partial night losses due to weather and instrument prob- 
lems, but this simple calculation serves our purpose). 
This value is consistent with the one cluster SN la spec- 
troscopically confirmed. It is difficult at this time to 
push further on this comparison because the cluster SN 
la (Abell 1246/ID10) was preferentially observed due to 
its proximity to the cluster center and clear offset from 
the host galaxy. 

A separate consistency check involves comparing the 
number of expected SN la with the number of IC SN can- 
didates (although it should be noted that all of our IC SN 
candidates are at R > r2oo while the number of expected 
cluster SN la was calculated for R < 7'2oo)- Assum- 
ing that all three of the IC SN candidates found in the 
cluster-centered chip 1 are actual IC SN la, then ~3/13 
would be the resulting hostless to total SN la fraction, 
implying a ^20% IC stellar mass fraction. This number 
is in broad agreem ent with direct measurements of the IC 
stellar population (Feldmeier et a l. 2004: Gonzale z et al.l 



I2005t IZibetti et al.ll2005t iKrick fc Bernstein! |2007[ ). We 
discuss the IC SN candidates further in § [9] 

8.2. Core Collapse SN 

Core collapse SN progenitors are massive, necessar- 
ily young stars, which is consistent with the negligibly 
small CC supernova rate in E/S0 type galaxies (e.g., 
iMannucci et al.ll2005l ). Therefore one might not expect 
CC SN to contribute significantly to our survey. In gen- 
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eral, CC SN are fainter than their SN la counterparts, 
leading to shorter visibility times. In detail this is not 
always the case since individual CC SN e can be quit e 
bright or have a very long plateau phase (|Hamuvll2003l) . 
which would prolong the visibility time. However, it so 
happens that, as we describe below, the CC SN rate in 
our survey is comparable to that of la's. 

We begin the exploration of this class by examin- 
ing results from previous surveys. The Mount Stromlo 
Abell cluster su pernova survey (| Germany et al.l 12004 
iReiss et aLlll998l ). which ran for 3.5 years, searched for 
SN in galaxy clusters at a redshift of 0.04 < z < 0.2. 
Of the 25 cluster SN discovered, 14 were SN la (identi- 
fied spectroscopically or by their light curve), while the 
remaining 11 are likely to be CC events (because the 
Mount Stromlo Abell cluster supernova survey obtained 
light curves for all their candidate SN, it is not likely 
that these events were AGN or varia ble stars). Alter- 
native ly, taking the 136 SNe sample of fCappcllaro et ahl 
(1999), and selecting the cluster SNe, iMannucci et al.l 
( 20071) found 44 cluster SNe, 16 of which are CC. These 
two results suggest that somewhere in the neighborhood 
of 40% of the detected cluster SNe may be CC. 

Empirically, no cluster SNe search has uncovered 
more CC SNe than SNe la (iGermanv et all 12004 
iGal-Yam etHI l2007t IMannucci et all 120071) . As such, 
we conservatively suggest that even if ^10 of our ex- 
cess cluster transients are due to CC SN, we still have an 
unexplained population of ^20-25 cluster transients. 

8.3. Cluster AGN 

Active galactic nuclei and quasars exhibit optical 
variability on a variety o f time scales, from days 
to years (iWebb fc Malkanl [20001 : iTotani et all [20051 : 
iKlesman fc Saraiedinil 120071 ). with somewhere between 
30-50% of all AGN displaying variability during the 
course of these surveys. Indeed, many of the transients 
found in this survey with known counterparts are QSOs 
(see § H]) . We now test whether the remaining detected 
central transients can be AGN. 

Multiwavelength studies of clusters find that roughly 
5% of cluster mem bers with Mr < —20 are AGN 
(jMartini et al.ll2007ft . These objects are more centrally 
concentrated than cluster g a laxies of the same luminos- 
ity, although iMartini et al.l (|2007f ) only probe the area 
covered by the ACIS-I chips on board Chandra, which 
is always R < r 2 oo for their samp le. Adoptin g a 5% 
AGN fraction and applying the iPopesso et alJ (|2007l ) 
A^oo vs M200 scaling relation (for M r < —20) we cal- 
culate that our sample contains ^200 cluster AGN. Of 
these, ~ 30 — 50% are variable implying that as many as 
60 to 100 cluster AGN could have been detected. 

There are various reasons why this number exceeds the 
actual number we are trying to account for, ~ 20 — 25. 
First, our temporal coverage is sparse and many of our 
clusters only have two imaging epochs so far. As the 
survey continues we expect to detect more of the cluster 
AGN population. Also, our estimate for the number of 
cluster AGNs relies heavily on galaxy cluster scaling re- 
lations based on different samples and techniques, which 
can plausibly introduce an error of a factor of ~2. Al- 
though this doesn't introduce a bias in one direction, 
it does result in making the difference between the pre- 
diction and observation less statistically significant. We 



conclude that it is plausible that the remainder of the 
statistical transient excess in clusters is due to AGN. In- 
deed, if we have over-estimated the number of cluster 
SNIa and CC SNe in the previous two sections (we have 
tried to get them correct at the factor of ^2 level), it is 
also plausible that cluster AGN can make up the differ- 
ence and account for the cluster transient excess. 

9. THE INTRACLUSTER SUPERNOVAE CANDIDATES 

As has been discussed, the number of IC SN candidates 
in the cluster-centered chip 1 is in rough agreement with 
the expected ~10 cluster SN la in the survey thus far 
(see § 15.11 and 18. ip . Even if the simplest explanation for 
the three IC SN candidates (all at R > ?"2oo) is that they 
are background events whose host was not detectable 
in our stacked images, it is interesting for the sake of 
argument to explore the implications of the three IC SN 
candidates being associated with the clusters. Note that 
the following discussion is not significantly altered if one 
takes the 'background' rate of hostless events to be that 
measured by chip 3, from which we would conclude that 
only two out of three of our IC SN candidates are actually 
associated with the cluster. 

A comparison between our IC SN candidates and ICL 
studies in the literature points to an untapped regime 
in IC studies. First the three IC SN candidates all lie at 
R > ^200 (or R > 1500 kpc) from the BCG. This is an un- 
explored radial regime among observational ICL studies, 
given that the study with the hi ghest radial measurement 
is that of IZibetti etaH (|2005h who stacked 683 SDSS 
clusters and was able to probe the aver age ICL out to 
~900 kpc. A simple reference point from Gonzale z et al.l 
(|200l can be used to show how unusual three IC SN at 
R > ^200 really is. From their sample of 24 clusters with 
direct low surface bri ghtness measurements of the ICL, 
iGonzalez et "all ((20071) find that 80% of the BCG+ICL 
light (of which the ICL dominates) is located within 300 
kpc of the cluster center, although extrapolation is nec- 
essary to extend this result out to R > r 2 oo- Using this 
as a baseline, finding three IC SN at R > r2oo would 
imply that we would expect at least 12 IC SN within 
the central 300 kpc of our clusters! The fact that none 
were found suggests that either the three IC SN are in 
fact not associated with the cluster, there is an excess 
of heretofore undetected ICL beyond r2oo, or that the 
stellar population in the cluster outskirts is more con- 
ducive to SN events. The answer could plausibly be 
a combination of the second and third possibilities just 
mentioned. For instance, infalling groups at the cluster 
outskirts would presumably have their own ICL (unac- 
counted for in studies at smaller radii), and it may be 
IC SN in these infalling groups that we are detecting. If 
these IC stars (whether or not they are associated with 
infalling groups) are more recently stripped than those in 
the central regions of the cluster (or if a sprinkle of star 
formation has been induced in the process of the stars 
becoming unbound), then these outskirt IC SN could be 
either CC events (which is compatible with their single 
epoch magnitude) or prompt SN la seen to be associated 
with star forming regions. In any case, the possibility of 
IC SN in the outskirts of clusters is extremely intrigu- 
ing and warrants a more coordinated effort. To investi- 
gate this potential large cluster radius IC SN population 
further requires continued cluster SN searches utilizing 
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wide-field imagers capable of imaging out to the cluster 
outskirts. 

It should also be noted that given that this initial 
phase of the survey expected to find ~10 cluster SN, 
some IC SN candidates within T2oo would be expected 
given that the ICL fraction typically seen in the liter- 
ature is ^20%. The fact that we found no candidates 
within r2oo is surprising, and a larger data set is needed 
to investigate if the true cluster SN la rate is lower than 
expected from previous cluster studies. For instance, 
if we take the early galaxy type SN la rate measure- 
ment from iMannucci et al.l (|2005h (with a SN la rate of 
0.038 SNuM) as our cluster SN la rate rather than of 
iSharon et al.l (|2007bf ) , then we would have expected only 
~4 cluster SN la in the survey thus far and plausibly 
no IC SN la. Further study is needed to resolve this 
situation. 

10. SUMMARY AND FUTURE WORK 

We have presented initial results from our program to 
search for and identify SNe la in a sample of X-ray se- 
lected clusters at 0.1 < z < 0.2. Of particular interest 
are hostless SN whose progenitors are IC stars. This 
survey will measure the SN la rate in galaxy clusters 
and the mean ICL fraction by noting the relative num- 
ber of hosted to hostless events. Ultimately, we aim to 
constrain the metal enrichment of the ICM and the role 
that IC SNe play. 

The data presented in this paper include our intitial 
photometry-only campaign followed by our first spectro- 
scopic results. We use the 30 x 30 arcminute field of view 
of Chip 1 on the 90Prime imager as our cluster-centered 
field, which allows us to probe beyond r2oo in ah of our 
clusters. We also use one of the off-cluster CCDs as a 
probe of the background rate of transients. 

The initial results presented in the paper include - 

1. Using our automated pipeline - which reduces, dif- 
ference images, and posts candidate transients to 
the web for human inspection - we have identified 
218 unique transient sources in the cluster-centered 
Chip 1 (107 of these were within ^200) and 102 on 
the off-cluster Chip 3. 

2. Of the transient sources we have identified, four 
were apparently hostless events (three of the four 
were in the cluster-centered chip 1). Continued 
imaging of these fields will put tighter constraints 
on possible faint hosts, while our current limits 
reach a magnitude limit of at least g ~25 (cor- 
responding to a galaxy of M g > —14.3 in all 
cases). Interestingly, the three IC candidates on 
the cluster-centered chip were all at R > r2oo, in- 
dicating that they are either not genuine cluster 
events or that some mechanism must be invoked to 
cause an excess of IC SN in the cluster outskirts. 
In any case, further study is needed. 

3. Initial spectroscopic results from one night of MMT 
Blue Channel Spectrograph data shows that we can 
identify and followup cluster SN la and other inter- 
esting variables. Of the four confirmed SNe la, one 
is in a cluster. Additionally, we tentatively identify 
a SNe Ib/c that could be in a cluster. Spectroscopic 



followup will become a major focus as the survey 
continues. 

4. Plotting the number of transients as a function of 
a scaled radius, -R/V200, reveals an excess of tran- 
sients within 1.25r2oo in our clusters. Of the ^40- 
50 excess transients within r2oo, we conservatively 
estimate that ~10 of these are cluster SN la, ~10 
are CC SN of some form, and the rest are cluster 
AGN. 

This survey is in its beginning stages and there are 
several avenues of future research. The main focus will 
continue to be on finding, and spectroscopically identi- 
fying, hostless and hosted cluster SN la. Inevitably, this 
will allow for other time domain and cluster science. For 
instance, our variability survey offers an opportunity to 
study cluster AGN out to a much larger radius than is 
feasible with most m ulti-object spectrogr aphs and X-ray 
instrumentation (e.g.. lMartini et aTl 2007). Using the full 
data set (all four CCDs) to search for all transients on 
long, ^monthly, time scales along with fast optical tran- 
sients (objects that vary on ^1000 sec time scales) and 
faint eclipsing systems are among our future plans. 

We are in the process of obtaining deep auxiliary r- 
band images of all of our cluster fields. The resulting 
red seqence identification will confidently identify old, 
early-type galaxies. Any transients associated with these 
galaxies are highly likely to be SN la, since CC SN rarely 
occur in old stellar populations. This will also be a 
unique, 0.1 < z < 0.2 data set for studying the cluster 
red sequence and the luminous red galaxy halo occupa- 
tion distribution. 

Measuring the SN rate in low redshift galaxy clusters 
(both from hosted and hostless systems) will provide a 
snapshot of the processes that contribute to the metal en- 
richment of the ICM. In combination with ongoing clus- 
ter SN research at higher redshift, we aim to constrain 
the SN la progenitor associated with old stellar popula- 
tions, and ultimately understand the metal enrichment 
history of the ICM. 
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Fig. 1. — Histogram of the distribution of stellar FWHM during our 90Prime imaging campaign. 
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Fig. 2. — Histogram of the distribution of image zeropoints with respect to the median of the distribution. 
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Fig. 3. — Transient point source detection efficiency as a function of input magnitude for three values of the seeing. The bright end 
cutoff in efficiency is due to our strict rejection of bright candidates. We are better able to detect bright transients during poorer seeing 
conditions simply because the transients are further from saturation. The detection efficiency has been rescaled to account for the area of 
the chip that was masked for any reason. 
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Fig. 4. — Transient point source detection efficiency as a function of zeropoint deviation from the median. To achieve high detection 
efficiency over the g ~19-21.5 region necessary for this survey, clear conditions are necessary. The detection efficiency has been rescalcd to 
account for the area of the chip that was masked for any reason. 
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Fig. 5. — The difference between the magnitude of the artificial transient placed into an image versus the recovered magnitude as a 
function of input magnitude. Shown are the results from a typical run as discussed in § [4] 
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Fig. 6. — Histogram showing the magnitude distribution of identified transients. The dashed histogram is plotted after the removal of 
transients that were detected as variables in multiple observing epochs. The dotted histogram is the transient magnitude distribution after 
removal of all sources known to not be at the cluster redshift. 
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Fig. 7. — Rest wavelength spectra of SNe la. Overplotted are SNe la template spectra from Peter Nugent (Nuecnt ct al3 l2002t ). 
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Fig. 8. — Spectra of six low redshift (2 < 0.3) galaxies. For a discussion, see § 16.41 
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Fig. 9. — Spectra of five high redshift QSOs. The location of the A-band atmospheric absorption feature is indicated with a telluric 
symbol. No attempt was made to remove this feature. 
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Fig. 10. — Spectra of five intermediate rcdshift QSOs. The location of the A-band atmospheric absorption feature is indicated 
telluric symbol. No attempt was made to remove this feature. 



with a 



Photometric Search for Transients in Galaxy Clusters 



23 



20 



15 



CO 

d 
o 
o 

+ 

0) 

o 

CO 



10 



5 







i i i i i i i i i i i i i i i i i i i i i i i i i r 




A1246/ID9 




A516/1D3^ 
"A620/ID5 




e 



.A620/ID12 




J I I I I I I I I I I I I I I I I I I I I I I I I L 



3000 4000 5000 6000 7000 8000 

Rest Wavelength (A) 



Fig. 11. — Spectra of five variable stars. The location of the A-band atmospheric absorption feature is indicated with a telluric symbol. 
No attempt was made to remove this feature. 
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Fig. 12. — Spectra of two unknown objects, although we suspect that the top spectrum may be a SNI b/c roughly eight days past explosion 
at a redshift of z = 0.1, with some obvious galaxy contamination. We overplot a template SN spectrum with these characteristics using 
the template of Nugent ct al. (2002). Even though Abell 516/ID10 was very faint at discovery, we devoted some spectroscopic resources 
to it since it was an apparent IC event. The location of the A-band atmospheric absorption feature is indicated with a telluric symbol. No 
attempt was made to remove this feature. 
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Fig. 13. — Radial profile of the variable sources as a function of clustcrccntric radius, scaled by r200- All values have been renormalized 
to account for the different areas probed by the radial bins (including masked regions). All of the clusters have been probed out to r200- 
The 'background' transient rate is represented by the horizontal bar and was measured from chip 3. 
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Fig. 14. — Examples of template g band SN la light curves at the redshift edges of our cluster sample (0.1 < z < 0.2). 
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TABLE 1 
List of Cluster Fields 



V^/l Lib LCi 


z cl 


us 




,5 




JZjJJUdl _L 














( T9finn c\\ 


( T9nnn c\\ 

^ J ZUUU.U J 


[±u ergs b 










A2703 


0. 


114 


00:05:30.6 


+16:5:11.9 


2.72 


11/18/06* 


12/15/06 






A136 


0. 


157 


01:04:14.4 


+25:05:03.9 


1.32 


11/17/06* 


12/15/06 






A351 


0. 


.105 


02:25:42.3 


-08:40:21.8 


0.21 


11/18/06* 


12/15/06 






Z0256 


0. 


.194 


02:59:32.6 


+00:13:56.1 


0.85 


11/17/06* 


12/15/06 






A516 


0. 
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04:50:24.9 


-08:49:24.6 


0.17 


11/18/06* 


12/15/06 
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0. 
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0. 
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0. 


124 


10:16:23.2 


+33:38:19.8 


3.12 


01/28/06 
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1.00 
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08:05:43.2 
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2.75 
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02/25/06 
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12/15/06 
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0. 
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+49:08:40.1 


6.71 


01/29/06 


02/25/06 


04/25/06 
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0. 
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11:27:01.8 


+17:07:47.7 


2.07 


01/29/06 


02/25/06 






A1246 


0. 
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11:23:58.8 


+21:28:49.0 


6.31 


01/29/06 


02/24/06 


12/15/06+ 




A1553 


0. 


.165 


12:30:48.9 


+10:32:48.4 


7.19 


01/29/06 


02/24/06 


04/23/06 




A 1204 


0. 
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11:13:20.5 


+17:35:41.0 


6.50 


01/29/06 


02/25/06 


04/24/06 




Z1432 


0. 
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07:51:25.0 
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5.18 


01/30/06 


04/25/06 
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0. 
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08:22:09.3 
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0. 
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08:25:28.8 
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0. 
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08:19:25.7 
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0. 
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09:24:05.4 


14:10:24.2 


5.70 
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02/26/06 






A923 


0. 
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10:06:38.9 


+25:54:44.5 


2.05 


01/30/06 


02/26/06 






A1033 


0. 


126 


10:31:44.3 


+35:02:28.3 


5.12 


01/30/06+ 


04/25/06+ 


12/15/06+ 




A1132 


0. 
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10:58:23.7 


+56:47:42.7 


6.76 


01/30/06 


04/25/06 






Aim 


0. 


165 


10:50:46.7 


-02:21:24.1 


4.12 


01/30/06 


02/25/06 






Z3179 


0. 


143 


10:25:58.1 


+12:41:10.6 


4.71 


01/30/06 


04/25/06 






RXJ1156 


0. 
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11:56:55.8 


+24:15:37.4 


3.10 


01/30/06 


04/25/06 




12/15/06+ 
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0. 
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11:22:16.0 


+19:35:49.3 


1.70 


01/30/06 


02/25/06 


04/25/06 


A1366 
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11:44:36.9 


+67:24:23.5 


3.91 


01/30/06 


02/25/06 


04/25/06 




A1902 


0. 


.160 


14:21:40.5 


+37:17:32.6 


4.55 


04/23/06 


06/19/06 






A1677 


0. 
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13:05:50.8 


+30:54:13.7 


5.37 


04/24/06 


06/19/06 
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0. 
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14:25:56.7 


+37:48:59.0 
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04/24/06 


06/19/06 
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0. 


131 


14:32:38.0 


+31:38:50.5 


3.99 


04/24/06 


06/19/06 
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0. 


.164 


17:20:09.9 


+27:40:08.9 


5.95 


04/25/06 


06/19/06 
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0. 
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16:35:49.6 


+66:12:44.8 


9.55 
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06/19/06 
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17:50:17.0 


+35:04:59.7 


5.34 


04/25/06 


06/19/06 







* Denotes dates when 90Prime readout noise was anomalously high. These epochs should not be con- 
sidered as part of the main sample, but are included for those clusters for which they are the only other 
imaging epoch. 

+ Chip 3 data for this epoch was not used either due to the presence of a vary bright field star in the field, 
poor focus, or an automated data reduction error 
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TABLE 2 

Summary of Intracluster SNe Candidates 



Cluster 


a 


S 


9 


Rbcg 


Rbcg/i~2oo 








(J2000.0) 


(J2000.0) 


[mag] 


[kpc] 








A516 


04:49:42.5 


-08:51:42.2 


23.83 ± 0.90 


1595 


2.0 


25.1 


-14.0 


A1201 


11:13:12.6 


+13:12:10.3 


22.52 ± 0.09 


2537 


1.4 


25.3 


-14.3 


A1413 


11:56:9.63 


+23:32:29.4 


22.12 ± 0.05 


2172 


1.0 


26.1 


-13.1 


A2259 


17:23:33.4 


+27:49:31.6 


22.54 ± 0.08 


7764 


4.2 


25.3 


-14.2 
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TABLE 3 

List of Identified Variables - Cluster Centered Chip 1 



Cluster 


Epoch 


a 


8 


9 


Dist from BCG 


Notes 






(J2000.0) 


J2000.0 


mag 


kpc 





A136 2006-12-16 02:43:01 

A136 2006-12-16 02:43:01 

A136 2006-12-16 02:43:01 

A351 2006-12-16 03:57:01 

A351 2006-12-16 03:57:01 

Z0256 2006-12-16 05:11:56 

Z0256 2006-12-16 05:11:56 

Z0256 2006-12-16 05:11:56 

Z0256 2006-12-16 05:11:56 

A516 2006-12-16 04:34:40 

A516 2006-12-16 04:34:40 

A516 2006-12-16 04:34:40 

A743 2006-04-25 05:12:57 

A743 2006-12-16 08:53:47 

A743 2006-12-16 08:53:47 

A743 2006-12-16 08:53:47 

A565 2006-01-28 02:42:49 

A565 2006-02-24 02:13:37 

A565 2006-04-23 03:04:32 

A565 2006-04-23 03:04:32 

A572 2006-02-24 02:48:27 

A572 2006-02-24 02:48:27 

A572 2006-02-24 02:48:27 

A572 2006-02-24 02:48:27 

A572 2006-02-24 02:48:27 

A572 2006-04-25 03:15:23 

A572 2006-04-25 03:15:23 

A580 2006-02-24 03:49:52 

A580 2006-04-24 03:05:59 

A580 2006-04-24 03:05:59 

A580 2006-04-24 03:05:59 

A580 2006-04-24 03:05:59 

A580 2006-04-24 03:05:59 

RXJ0759 2006-01-28 05:07:38 

RXJ0759 2006-01-28 05:07:38 

RXJ0759 2006-01-28 05:07:38 

RXJ0759 2006-04-24 03:52:36 

RXJ0759 2006-04-24 03:52:36 

RXJ0759 2006-12-16 09:36:16 

RXJ0759 2006-12-16 09:36:16 

RXJ0759 2006-12-16 09:36:16 

RXJ0759 2006-12-16 09:36:16 

RXJ0759 2006-12-16 09:36:16 

RXJ0804 2006-01-28 06:15:30 

RXJ0804 2006-02-24 05:20:32 

RXJ0804 2006-02-24 05:20:32 

RXJ0804 2006-02-24 05:20:32 

RXJ0804 2006-02-24 05:20:32 

RXJ0804 2006-04-24 04:33:50 

RXJ0804 2006-12-16 07:07:14 

RXJ0804 2006-12-16 07:07:14 

RXJ0804 2006-12-16 07:07:14 

RXJ0804 2006-12-16 07:07:14 

RXJ0804 2006-12-16 07:07:14 

RXJ0804 2006-12-16 07:07:14 

A667 2006-01-28 06:52:55 

A667 2006-01-28 06:52:55 

A667 2006-01-28 06:52:55 

A667 2006-02-24 05:57:42 

A667 2006-04-24 05:11:06 

A667 2006-04-24 05:11:06 

A750 2006-04-24 05:49:13 

A853 2006-01-28 08:36:18 

A853 2006-01-28 08:36:18 

A961 2006-02-24 07:52:05 

A1201 2006-01-28 10:18:32 

A1201 2006-01-28 10:18:32 

A1201 2006-04-24 07:09:10 

A1201 2006-04-24 07:09:10 

A1201 2006-04-24 07:09:10 

A1413 2006-01-28 11:09:51 

A1413 2006-01-28 11:09:51 

A1413 2006-02-24 09:06:07 

A1413 2006-02-24 09:06:07 

A1413 2006-02-24 09:06:07 

A1413 2006-04-24 07:59:00 

A1413 2006-04-24 07:59:00 

A1413 2006-04-24 07:59:00 

AlfiSQ 9n0fi_ni_98 11-M'IQ 



1:3:19.9 
1:4:2.14 
1:3:55.5 
2:26:21.3 
2:25:56.9 
2:59:35.2 
2:59:7.80 
3:0:18.8 
3:0:15.3 
4:49:42.5 
4:50:52.5 
4:50:29.1 
9:7:25.7 
9:6:55.2 
9:6:32.4 
9:7:12.5 
7:7:26.8 
7:9:33.9 
7:10:10.5 
7:7:49.4 
7:15:26.0 
7:14:11.5 
7:13:8.08 
7:12:56.3 
7:12:38.0 
7:13:46.5 
7:13:25.1 
7:26:15.1 
7:25:5.34 
7:26:22.4 
7:26:20.5 
7:26:15.1 
7:26:12.5 
8:0:37.4 
7:59:17.3 
7:58:34.4 
7:59:40.3 
7:58:32.1 
8:0:46.0 



52.3 
13.0 
45.9 
32.1 



7 
7 
7 

7: 

8:5:10.2 
8:5:6.84 
8:4:28.0 
8:4:24.6 
8:4:8.89 
8:2:51.8 
8:5:18.5 
8:3:24.0 
8:3:12.9 
8:4:50.2 
8:4:21.4 
8:3:57.0 
;:28:49.1 
■:28:31.9 
■:28:15.9 
:27:9.95 
:28:49.1 
1:26:41.5 
9:9:47.8 
9:43:12.5 
9:42:53.6 
10:17:29.6 
11:13:36.0 
11:12:17.8 
11:13:29.0 
11:13:12.6 
11:12:47.5 
11:55:50.3 
11:55:4.32 
11:56:9.63 
11:56:5.11 
11:55:46.6 
11:55:55.2 
11:55:44.4 
11:55:42.0 
1 a.in-Rn 9 



25:4:28.4 
25:7:22.8 
24:50:20.5 
-8:54:19.8 
-8:45:38.3 



0:18 
0:22 
0:19 
0:13 

-8: 



:38.6 
:33.9 
:11.6 
;:28.4 
1:51:42.2 
-8:45:3.96 
-8:50:2.11 
10:24:30.9 
10:14:13.6 
10:9:29.1 
10:23:34.9 
71:44:42.2 
71:34:49.1 
71:50:20.4 
71:52:31.6 
54:37:29.0 
54:26:52.5 
54:26:39.5 
54:28:31.2 
54:41:0.311 
54:38:22.6 
54:42:15.0 
41:37:34.2 
41:32:48.4 
41:24:43.2 
41:14:48.2 
41:37:34.2 
41:22:34.1 
53:47:50.4 
54:2:32.0 
54:10:30.8 
53:55:49.2 
54:6:42.1 
53:47:15.5 
53:47:44.5 
54:4:42.2 
54:1:17.4 
54:6:42.9 
46:43:32.8 
46:35:19.0 
46:48:15.8 
46:56:55.1 
46:39:10.9 
46:48:42.9 
46:35:40.3 
46:55:18.9 
46:58:12.1 
47:1:45.2 
46:57:15.4 
46:40:37.7 
44:44:38.0 
44:34:3.13 
44:47:24.2 
44:29:25.1 
44:44:38.0 
44:55:37.4 
11:6:36.6 
15:35:9.70 
15:26:4.16 
33:45:15.6 
13:38:31.0 
13:21:6.03 
13:34:42.0 
13:12:10.3 
13:14:4.21 
23:23:58.5 
23:31:18.7 
23:32:29.4 
23:11:16.5 
23:24:47.8 
23:15:44.3 
23:37:29.5 
23:28:55.6 

_1 -OA-OD 7 



24.51 ± 0.60 
22.91 ± 0.12 
21.75 ± 0.06 
22.71 ± 0.17 
21.71 ± 0.07 
22.87 ± 0.14 
22.81 ± 0.14 
23.59 ± 0.31 
22.57 ± 0.14 
23.83 ± 0.90 
22.09 ± 0.07 
22.93 ± 0.19 
21.54 ± 0.06 
22.00 ± 0.09 
21.32 ± 0.05 
22.29 ± 0.12 
22.14 ± 0.09 

21.52 ± 0.05 
21.56 ± 0.07 
21.93 ± 0.09 



21.80 
23.12 
22.93 
21.03 
20.90 
22.15 
23.23 
21.71 



0.04 
0.15 
0.14 
0.03 
0.02 
0.06 
0.18 
0.04 



22.15 ± 0.07 
21.87 ± 0.06 
20.40 ± 0.02 
20.76 ± 0.02 
22.20 ± 0.08 
22.01 ± 0.06 



22.19 
21.93 
22.64 
22.71 
22.59 
21.27 
22.20 
22.09 
22.71 
21.98 
22.58 



0.07 
0.06 
0.10 
0.10 
0.07 
0.02 
0.04 
0.04 
0.06 
0.04 
0.07 



22.40 ± 0.06 
22.08 ± 0.04 
19.13 ± 0.00 
23.72 ± 0.32 
21.60 ± 0.03 
21.95 ± 0.06 
21.90 ± 0.05 
21.70 ± 0.05 



22.83 
22.32 
20.47 
22.24 
22.29 
24.18 
20.73 
23.25 
22.34 
22.00 



0.11 
0.07 
0.01 
0.06 
0.07 
0.47 
0.02 
0.22 
0.09 
0.13 



22.00 ± 0.12 
22.83 ± 0.12 
22.57 ± 0.10 

21.55 ± 0.04 
22.87 ± 0.12 
22.52 ± 0.09 
21.95 ± 0.05 
22.49 ± 0.08 
22.19 ± 0.07 
22.12 ± 0.05 
23.64 ± 0.24 
21.59 ± 0.04 

21.56 ± 0.03 
22.29 ± 0.07 
22.03 ± 0.06 
i s -4- n on 



2012 
589 
2498 
1961 
740 
919 
2054 
2455 
2069 
1595 
1206 
183 
1620 
775 
1425 
1158 
528 
1622 
2077 
1429 
1188 
1628 
2025 
1984 
1666 
615 
897 
1897 
1753 
640 
1204 
1897 
375 
1728 
466 
1607 
528 
1362 
1866 
1466 
680 
934 
1363 
1716 
2672 
292 
1825 
1561 
2890 
2845 
2386 
3010 
2881 
1885 
1465 
1177 
1958 
329 
2966 
1177 
2716 
1991 
3140 
1672 
2063 
2759 
1779 
2070 
2537 
2114 
1112 
1173 
2171 
2537 
989 
1808 
2198 
1088 



z=0.172 galaxy; BCS 

z=0.246 galaxy; BCS 
SN 2006ub; SNIa; BCS 
SN 2006ua; SNIa; BCS 

z=0.199 CG; BCS 



?; BCS; IC Cand 
Stellar spec; BCS 
Stellar spec; BCS 

z=0.267 Galaxy; BCS 
SN 2006ud; SNIa; BCS 



RD 

z=1.803 QSO; SDSS 
RD 

CG; SDSS 



RD 



z=2.93 QSO; BCS 
SN 2006uc; SNIa; BCS 
RD 



z=1.775 (phot) QSO/1 
z=1.775 (phot) QSO/1 

RD 

RD 



z=0.231 QSO; SDSS 
IC Cand 

QSO z=1.524; SDSS 



IC Cand 



tin riTTTTST nssnAQ/*? 



Sand et al. 



TABLE 4 

List of Identified Variables - Background Chip 3 



Cluster 


Epoch 


a 

(J2000.0) 


8 

J2000.0 


g 
mag 




Dist from BCG 
kpc 


Notes 


R Y T1 7^0 

JTwV.J 1 I OU 


2006-04-26 


11:26:39 




OO.O.O 1 . 1 


20.95 


± 


n 04 


oyzo 




pY T1 7^0 


2006-04-26 


11:26:39 


1 7- r >9-4fi ^ 

J. t . O £ .1U .O 


oo. i . y . i u 


21.72 


± 


0.08 


5349 




A2218 


2006-04-26 


10:46:38 




fifi-99-48 9 


21.35 


± 


0.03 


5670 




A2218 


2006-04-26 


10:46:38 


1 fi-40-49 & 


fifi-1 1 -10 Q 


22.32 


± 


0.07 


5302 




A2259 


2006-04-26 


10:10:35 


17:23:33.4 


97-4Q-31 fi 


22.54 


± 


0.08 


7764 


TP PanH 


A2259 


2006-04-26 


10:10:35 


1 7-99-7 9D 


97-9fi-9Q 7 


21.87 


± 


0.10 


4963 




A1677 


2006-06-19 


04:59:20 


13:9:27.7 


30-44-19 4 

OU . 11 . li/trr 


20.54 


± 


0.02 


8794 




A1677 


2006-06-19 


04:59:20 


13:9:3.27 


31 -fi-9fi fi 
Ol.U. zu.u 


22.18 


± 


0.07 


7930 




A1689 


2006-01-28 


11:51:19 


1 3-1 4-46 1 


-1:8:52.2 


21.92 


± 


0.06 


9305 




A1689 


2006-01-28 


11:51:19 


13:13:46.3 


-1:13:19.2 


22.07 


± 


0.06 


6431 




A1689 


2006-02-24 


09:46:20 


-LO. -LO .01.0 


i -i r-^ no 

- J. . ±0.0 . UU 


22.17 


± 


0.07 


5781 


Zj — J. .U 1 O ^ jJLLOLr J Vc^OV_7 j -L 


A1689 


2006-06-19 


04:03:19 


1 3-1 4-^0 8 

_LO . 1t:.iJU.U 


-1 -29-45 9 

_L . £i C7 . rtO . v? 


20.82 


± 


0.03 


8514 


A1689 


2006-06-19 


04:03:19 


13:14:17.9 


-1:35:12.1 


22.13 


± 


0.10 


8208 




A1689 


2006-06-19 


04:03:19 


1 3-1 3-48 7 

-LO . -LO . lO . t 


-1 ^^i-^Q ^ 

- ± . £0 'OlJ >o 


22.12 


± 


0.09 


6472 




A1902 


2006-04-24 


10:42:43 


1 4-94-37 9 

_L4:.Z4:.0 I . Z 


O / .ZU.OO.'i 


22.76 


± 


1 fi 

U. -LU 


^S^fi 
OOOU 




A1902 


2006-04-24 


10:42:43 


1 4-9fi-9 1 fi 

_L1. ZU .Z. 1U 


07.0.4 

O ( . O.rt. OO 


21.92 


± 


0.07 


8965 




A1914 


2006-04-25 


10:31:46 


1 4-9Q-97 7 


^7-4^-91 

O 1 . IO . Z _L . U 


21.98 


± 


0.05 


7316 


zj — j. .uy v^ou, O-L/oo 


A1914 


2006-04-25 


10:31:46 


1 4-98-99 1 

.L1.ZO . . _L 


1 .00 . -Li.o 


22.96 


± 


0.10 


5264 




A1914 


2006-06-19 


06:24:29 


1 4-30-Q fi3 


^8-0-41 Q 

OO .U.rt-L . y 


22.19 


± 


0.05 


8942 




A 1 cno 
-riiyou 


2006-06-19 


07:06:52 


1 4-3^-0 ^0 

j.'i.oo . y . ou 


Ol .O-L .'iO.O 


21.86 


± 


04 


4624 


Zj — U.OU Vo^DW, O-L/OO 


A 1930 


2006-06-19 


07:06:52 


14:35:0.177 


^1 -41 7 

OL.4L. OO. 1 


22.97 


± 


0.17 


4250 




A1930 


2006-06-19 


07:06:52 


1 4-34- r » c l Q 


^1 ^Q-94 ^ 
± . oy . zl 


21.93 


± 


0.05 


4108 




A1553 


2006-04-24 


10:03:28 


1 9-34-1 9 ^ 

_L Z . Ol . _L Z . O 


-Lu . zu . zo . y 


22.71 


± 


0.16 


8565 




J. lAtl Ki'iJO 


2006-01-29 


05:54:09 


10:1:32.4 


23-40-0 701 

iiO.rlU.U. 1 U-L 


22.28 


± 


0.10 


7137 




RY T0958 


2006-01-29 


05:54:09 


10-1 -25 3 

J_ O • J. . iiO . O 


23-56-20 2 

J . OU . ^U . 


21.05 


± 


0.03 


6955 




RX T0958 


2006-01-29 


05:54:09 


10:0:49.4 


93-54-33 7 

iiO . Ort. OO . 1 


21.82 


± 


0.06 


5664 


RD 


jwv.j uyuo 


2006-01-29 


05:54:09 


10:0:7.77 


9Q-44-c;fi 
ZO . H. OU . u 


22.60 


± 


0.12 


4078 




RY T0Q^8 
jwv.j uyuo 


2006-02-25 


05:03:22 


1 0-0-4Q 4 

J.U. U. iy . 1 


ZO . iJI.OO. 1 


21.62 


± 


0.05 


5664 


RD 


A1361 


2006-01-28 


12:29:23 


11:47:27.2 


4fi-94-^7 R 
1U . ZL / . u 


22.30 


± 


0.06 


4998 




A1361 


2006-02-24 


10:30:22 


1 1 -48-37 

_L J. . lO .O i . U 


4fv^4-1fi 

1U . OL -LU.O 


22.58 


± 


0.08 


6694 




A1361 


2006-02-24 


10:30:22 


11:47:7.16 


4fi-^4-9Q 9 

t:U . Ort. Zj . Z 


22.24 


± 


0.05 


4825 




A1361 


2006-04-24 


09:20:33 


11:47:7.34 


46:7:12.0 


23.04 


± 


0.11 


4909 




A1361 


2006-04-24 


09:20:33 


1 1 -47-4 fi 1 ^ 

LL.1I .1. UO 


4R-91 -49 


22.42 


± 


0.05 


4490 




A1413 


2006-01-28 


11:09:51 


1 1 ■ ^8-1 ^ 1 

_L 1 . Jo .LJ.l 


OQ.i /j. -49 1 

ZO . _L i ±.'±Z . 1 


22.74 


± 


12 


fi9Q1 
uzy 1 




A1413 


2006-04-24 


07:59:00 


1 1 ■ ^8-1 Q 8 

_L i . jo.iy.o 


OQ.Q1 • Q Q 1 
Z.O.OI-.OO. 1 


21.48 


± 


03 

u.uo 


fi^74 
uo I *± 




A1413 


2006-04-24 


07:59:00 


1 1 -^7-48 Q 

1L.JI .10 . U 


9^-1 ^-47 ^ 

ZO . -LO .11 .U 


21.12 


± 


0.02 


5450 




A961 


2006-02-24 


07:52:05 


1 0-1 8-^8 7 

_LU . -LO .OO. f 


^■47-1 ^ 8 

OO .11 .±0.0 


22.41 


± 


0.08 


4477 




A961 


2006-02-24 


07:52:05 


1 0-1 8-^1 9 


OO . IO .Ol.U 


21.05 


± 


0.02 


4222 






2006-04-24 


06:26:29 


1 n. 1 Q. c: cr o 
_LU. ly.oo.o 


Qq./ir-.c OO 
00.4:0.0.00 


22.02 


± 


0^ 
u.uo 


oy f 1 




A961 


2006-04-24 


06:26:29 


1 0-1 Q-48 4 
lu. 


^^-9fi-9 OQ 
00 .zu . z .uy 


22.98 


± 


0.12 


5948 




A961 


2006-04-24 


06:26:29 


10:19:23.1 


^■^9-1 4 

OO . OZ . _LU . 1 


22.59 


± 


0.08 


5075 




A853 


2006-01-28 


08:36:18 


9:45:19.1 


15:14:37.7 


21.71 


± 


0.09 


7665 




A853 


2006-01-28 


08:36:18 


Q-4^-4 3Q 


1 ^ - 1 9-47 Q 

_LO . -LZ .1 1 . iy 


22.22 


± 


0.14 


7143 




A853 


2006-01-28 


08:36:18 




1 ^-1 9-^7 1 

_LO . -LZ . O 1 . -L 


19.44 


± 


0.01 


6611 




A853 


2006-01-28 


08:36:18 


Q-44-3fi 3 

i7 .11. OU .O 


1 5-22-24 2 


22.29 


± 


0.29 


5782 




A853 


2006-01-28 


08:36:18 


9:44:10.3 


15-23-54 8 


21.17 


± 


0.06 


4720 




A853 


2006-04-25 


05:53:14 


9-44-35 5 

i7 .11. OO .O 


1 5-9-31 2 


21.93 


± 


0.08 


6186 




A667 


2006-04-25 


05:53:14 


8-31 -29 9 

o .o J. .iiu . <y 


44-39-30 7 

H. Oil .OU . 1 


22.24 


± 


0.06 


5904 


RD 


A667 


2006-04-24 


05:11:06 


8-32-24 1 

O . Oil . ill . _L 


44-53-53 fi 

n. 00 .00 . u 


22.66 


± 


0.11 


7063 


Zi — j-.OiiO 1 ^jiluii j Vc^kjv_// j. 


A667 


2006-04-24 


05:11:06 


8-31 -29 9 

o .o J. . . <y 


44-39-30 7 

rtrt.Oi .OU . 1 


22.25 


± 


0.08 


5905 


RD 


A580 


2006-01-28 


04:25:28 


7-29-51 

1 .£jZ> .O-L .U 


41 -32-55 8 


21.38 


± 


0.03 


5735 


RD 


A580 


2006-01-28 


04:25:28 


7-28-53 6 

I .iiO . OO .u 


41 -34-3fi 7 

1 J. . OLOU . 1 


22.52 


± 


0.07 


4470 




A580 


2006-04-24 


03:05:59 


7:29:51.1 


41 -32-56 3 

1 J. . Oil .OU . 


19.74 


± 


0.01 


5737 


RD 


A580 


2006-04-24 


03:05:59 


7:29:17.8 


41 -^8-9 40 
11.00 .Z.1U 


21.29 


± 


0.03 


5161 


LM V OO SU L1L / Zi 


A580 


2006-04-24 


03:05:59 


7-98-4fi Q 
i .zo . iu . y 


1 J. . Z c7 . Ol . U 


21.52 


± 


0.04 


4156 


A646 


2006-01-30 


02:51:19 


8:25:49.5 


47-1 0-97 ^ 

11 . _LU.Z / .O 


18.60 


± 


0.01 


5210 


Zi — U. -LO V_.' \_T iVVJl < / O 


A646 


2006-12-16 


10:22:04 


O . ZrO . -LU.O 


47-1 fi-^4 ^ 

1 1 . .LU .Ol. O 


21.86 


± 


0.07 


4498 


-7 — 1 49^ frtVint^ O^IO/I 

Zj — J..1ZO ^jJLLOL^ Vo^OV_/^ J. 


A 743 


2006-12-16 


10:22:04 


9:10:13.1 


1 0-4-9^ ^ 

.LU. 1. ZO.O 


21.97 


± 


0.12 


7557 


A 743 


2006-12-16 


10:22:04 


9:9:27.2 


1 0-^-8 98 


21.67 


± 


0.07 


5998 




A 743 


2006-12-16 


08:53:47 


9:8:58.4 


1 0-9^-^9 


21.05 


± 


0.04 


4743 




R Y T1 000 

llAtl ± OUU 


2006-02-25 


05:39:58 


1 0-^-7 00 

J-O. O . 1 . uu 


44-0-39 8 


21.08 


± 


0.03 


8018 




rx n nnn 


2006-04-25 


07:19:59 


10-5-3 80 


44:17:23.7 


21.85 


± 


0.06 


7906 




R Y T1 000 

JX-V-J 1UUU 


2006-04-25 


07:19:59 


10:3:17.2 


44:1:52.6 


21.75 


± 


0.05 


4890 




R Y T0804 

± UW WOU1 


2006-01-28 


06:15:30 


8:9:14.8 


46-39-38 9 

1U . OC7 .OO . <J 


20.86 


± 


0.02 


9560 


z =0.975 (phot) QSO/1 


RX T0804 

± 1.AJ UUU4 


2006-04-24 


04:33:50 


8:7:44.3 


46-38-44 

iu . 00 .11. vj 


20.56 


± 


0.02 


6729 


RX T0804 

-1 lAd UOyi 


2006-04-24 


04:33:50 


8:7:18.7 


47:1:18.1 


22.82 


± 


0.11 


6259 




A750 


2006-01-28 


07:59:36 


9:11:3.81 


11:4:28.3 


21.65 


± 


0.05 


4689 


7 _n o9« OSO- SDSS 

Zi U.fiiO Vi^OV./} O.L-'OO 


A750 


2006-02-24 


06:35:19 


9:11:56.1 


11:9:16.5 


22.53 


± 


0.07 


6974 




A750 


2006-04-24 


05:49:13 


9:11:36.7 


11:11:11.7 


22.71 


± 


0.12 


6308 




A750 


2006-04-24 


05:49:13 


9:11:33.0 


10:57:3.77 


21.90 


± 


0.05 


5794 




A750 


2006-04-24 


05:49:13 


9:11:22.5 


11:8:58.3 


22.91 


± 


0.14 


5632 




A750 


2006-04-24 


05:49:13 


9:11:5.82 


10:58:12.3 


22.24 


± 


0.07 


4665 




RXJ0741 


2006-04-25 


03:53:15 


7:44:23.5 


26:8:3.07 


21.35 


± 


0.05 


6684 


RD 


RXJ0741 


2006-01-29 


02:11:06 


7:43:56.4 


26:9:21.5 


21.94 


± 


0.06 


5784 




RXJ0741 


2006-01-29 


02:11:06 


7:44:46.9 


25:48:10.3 


20.83 


± 


0.02 


7574 


z=0.141 CG QSO; SDSS/4 


RXJ0741 


2006-02-25 


02:26:30 


7:44:49.8 


25:46:11.8 


20.87 


± 


0.02 


7756 


RXJ0741 


2006-02-25 


02:26:30 


7:44:23.5 


26:8:3.07 


22.22 


± 


0.08 


6685 


RD 


A620 


2006-01-29 


02:55:40 


8:9:43.3 


45:32:35.2 


20.21 


± 


0.01 


6157 





AK9n 9nnfi_n9_9c; m-n^ AQ s-s-9n 9 k-w-w n 99 nn -t- n or *?Q7i 



Photometric Search for Transients in Galaxy Clusters 



TABLE 5 

Blue Channel Followup of Transients 





Cluster/ID 


a 




5 


UT Date 


Classification 


z 


Discovery 


Exp. 




(J2000.0) 




J2000.0 


2007 Dec 20 






(mag) 


(sec) 


Ahpii 1033/TD0 


1 0-S2-40 


6 


4-S4-58-35 8 


0.54 


BAT, OSO 


1.56 


09 64+0 1 


600 


Abcll 1033/ID6 


10:30:31 


9 


+34:56:47.0 


0.53 


SNIb/c? 


0.1? 


23.16+0.19 


600 


Abcll 1246/ID2 


11:24:49 


9 


+21 -36 05 2 


0.51 


QSO 


0.92 


21.45+0.02 


600 


Abell 1246/ID9 


11:24:02 


5 


+21:36:58.1 


0.52 


Stellar spec 


0.0 


22.54+0.06 


600 


Abell 1246/ID10 


11:23:49 


6 


+21:28:40.4 


0.50 


Cluster SNIa 


0.19 


20.63+0.01 


900 


Abcll 136/ID3 


01:03:55 


5 


+24:50:20.6 


0.08 


Galaxy 


0.246 


21.75+0.06 


1800 


Abell 136/ID7 


01:03:19 


9 


+25:04:28.4 


0.10 


Galaxy 


0.172 


24.5+0.6 


1800 


Abcll 351/ID1 


02:26:21 


3 


-08:54:19.8 


0.15 


SNIa 


0.33 


22.71+0.17 


1800 


Abell 351/ID5 


02:25:57 





-08:45:38.4 


0.13 


SNIa 


0.27 


21.71+0.07 


900 


Abell 516/ID3 


04:50:52 


5 


-08:45:03.5 


0.23 


Stellar spec 


0.0 


22.09+0.07 


1800 


Abell 516/ID8 


04:50:29 


2 


-08:50:01.7 


0.26 


Stellar spec 


0.0 


22.93+0.19 


1800 


Abcll 516/ID11 


04:49:42 


6 


-08:51:43.1 


0.21 


?; IC Candidate 


? 


23.83+0.90 


1800 


Abcll 620/ID5 


08-06-14 


8 


+45-35-45 5 


0.29 


Stol 1 $\ V CiDPP 


0.0 


21 19+0 04 


600 


Abcll 620/ID12 


08:05:06 


9 


+45:32:37.3 


0.30 


Stellar spec 


0.0 


21.70+0.06 


1800 


Abell 646/ID13 


08:22:07 





+46:56:00.5 


0.32 


QSO 


2.36 


21.03+0.07 


1800 


Abell 646/ID19 


08:21:05 


6 


+47:16:33.0 


0.35 


QSO 


2.22 


21.40+0.10 


900 


RXJ0759/ID2 


07:59:13 





+54:04:42.2 


0.36 


QSO 


2.93 


22.20+0.04 


900 


RXJ0759/ID3 


07:58:45 


9 


+54:01:17.5 


0.38 


Galaxy 


0.21 


22.09+0.04 


1350 


RXJ0821/ID4 


08:21:23 





+07:39:41.8 


0.45 


QSO 


0.90 


22.20+0.11 


900 


RXJ0821/ID8 


08:20:52 


8 


+07:52:38.4 


0.40 


QSO 


0.795 


21.96+0.08 


900 


RXJ0821/ID9 


08:20:44 


6 


+07:50:07.6 


0.43 


QSO 


2.17 


22.43+0.13 


1350 


RXJ0821/ID15 


08:20:34 


3 


+07:39:38.1 


0.42 


Galaxy 


0.22 


21.73+0.07 


1350 


Z0256/ID19 


02:59:07 


8 


+00:22:33.0 


0.17 


Cluster Galaxy 


0.199 


22.81+0.14 


1800 


Abell 743/ID10 


09:06:55 


3 


+10:14:13.7 


0.48 


AGN 


0.267 


22.00+0.09 


900 


Abell 743/ID12 


09:06:32 


5 


+10:09:28.7 


0.46 


SNIa 


0.175 


21.32+0.05 


1350 



